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Abstract 

The self-consistent separable random-phase approximation (SRPA) model with Skyrme forces is extended 
to the case of magnetic excitations and applied to the description of spin-flip and orbital Ml giant resonances 
in the isotopic chain 142 ~ 152 Nd. The Skyrme forces SkT6, SkM*, SLy6 and SkI3 are used. The calculations 
show the onset of the scissors mode with increasing deformation. A specific three-peak structure of the 
spin-flip response is found and explained by particular neutron and proton spin-flip transitions. Although 
the employed forces provide an acceptable qualitative description, the Skyrme functional still needs further 



1 — 1 

improvement to reproduce quantitatively the experimental data. 



I. INTRODUCTION 



Magnetic giant resonances (GR) represent an important part of the nuclear dynamical response 
They were widely investigated within various phenomenological and microscopic models, for 
surveys see . Now more elaborated Skyrme, Gogny, and relativistic mean-field approaches 

^ u] based on density functional theory (DFT) are at our disposal. They are actively used for 
the description of nuclear ground states and dynamics covering mainly electrical modes but not 
magnetic ones. However magnetic modes could be very useful for analyzing and improving these 
approaches in the parts related to spin densities. Besides, since these modes (especially the spin- 
flip Ml GR) are sensitive to the spin-orbit splitting, they can be used to study the spin-orbit 
interaction in Skyrme and Gogny forces. They can also help to clarify the role of tensor forces and 
related spin-orbit densities [all], [lj| • 

This paper is devoted to an investi gati on of Ml spin-flip and orbital GR within the Skyrme- 
Hartree-Fock (SHF) approach [ill . I12I Earlier SHF studies of these GR were limited to a 



few explorations 



14 



151 ] and even these were not fully consistent. In the study 



a 



a hybrid 



model with partial inclusion of SHF in the Landau-Migdal formulation was exploited while the 
work [15|] used early Skyrme forces and omitted the important spin density. Quite recently a 



ft rst fully self-consistent systematic SHF investigation of the spin-flipMlGR was performed ^ 
within the separable Random-Phase- Approximation (SRPA) model [la, LL7(] extended to magnetic 



excitations 



id, 



la ]. It was shown that none of 8 different Skyrme parameterizations was able to 



describe simultaneously the one-peak structure in doubly-magic nuclei together with the two-peaks 
in deformed nuclei. 

In this paper, the SRPA for magnetic excitations is presented in more detail and applied to an 
investigation of spin-flip Ml GR in the chain of isotopes 142_152 Nd. The trends with increasing 
deformation and the number of neutrons are explored. Besides, we present our first results for 
the orbital scissors mode [3]. The calculations employ the Skyrme parameterizations SkT6 [3], 



SkM* 



ZM, SLy6 



2g| , and SkI3 



23]. 



II. SRPA FOR MAGNETIC EXCITATIONS 



General formalism 



SRPA is a fully self-consistent DFT model where both the static mean field and residual in- 
teraction are derived from the same functional. The present nuclear application is based on the 



2 



Skyrme functional 



11 



and then ex- 



tended to magnetic modes 



12l . 13 . It was first derived for electric excitations 

, [3] . The self-consistent factorization of the residual interaction in 
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SRPA considerably reduces the computational expense while maintaining a high accuracy. This 
makes the model very suitable for systematic studies. The SRPA residual interaction includes all 
contributions from the initial Skyrme functional as well as the Coulomb (direct and exchange) and 
pairing (at BCS level) terms. The model was widely used for the investigation of electrical GR in 
spherical and deformed (heavy and super-heavy) nuclei 



Starting point is the Skyrme functional with the energy density 
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where a, 71, bi, b[ are the force parameters. This functional involves time-even (nucleon p q , 

kinetic-energy r q , spin-orbit J q ) and time-odd (current j^, spin s q , and vector kinetic-energy T q ) 
densities where q denotes protons and neutrons. Densities without index, like p = p p + p n , are 
total. The contributions with bi (i=0,l,2,3,4) and b[ (i=0,l,2,3) are the standard terms responsible 
for ground state properties and electric excitations of even-even nuclei . In the standard SHF, 
the isovector spin-orbit interaction is linked to the isoscalar one by 6 4 = b^. The tensor spin-orbit 
terms tx 7t are often skipped. They can be switched in ([I]) by the parameter 71- • The spin terms 
with bi,b[ are relevant for odd nuclei and magnetic modes in even-even nuclei. Though b^b\ may 
be nni q nely deternnned as function, of M fl, then values were not y et weU tested by nuciear 
data. As was shown N.jnst these spin tonus are of a_nt importance fo r the spin-hip ML 
The general SRPA formalism is given elsewhere [16|, 113] • We present here only the basics and 
peculiarities of magnetic excitations. The SRPA simplifies the residual interaction of Skyrme RPA 
in a factorized (separable) form as 

1 K 

Kes P = ~q {Kqk,q'k'XgkX q 'ki + T] q k,q'k>Y q kY q >k'} (2) 

qq' k,k' = l 

where the indices q and q' label neutrons and protons, k numbers the separable terms, n q k,q'k' 



and r/gfc,g'A;' are the corresponding strength matrices, and and Yqk are time-even and time-odd 
hermitian one-body operators. We need these two kinds of the operators since the relevant Skyrme 
functionals involve both time-even and time-odd densities, see 



1(1 Q,Q. 



The model uses the Skyrme functional E[Jq(r,t)] = J drTisk( r ,t) depending on the local den- 
sities Jg = (p g , T q , J q , j g , Sg, T g) . The separable operators and strength matrices in (J2J) are self- 
consistently derived from this functional and read jlijj [l?}] 
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Here are the operators associated with the local densities Jq. Further, Q q k(r) and P q k(r) are 
generalized coordinate (time-even) and momentum (time-odd) hermitian one-body input operators 



which serve as generators of the separable terms 



16, 



171 ]. For EA modes, the input operators Q g k(r) 



are chosen and P q k{?) = i[H, Qqk] where H stands for the full Hamiltonian. Instead, for MA modes, 
the input operators P q k(r) are chosen and Q q k(r) = i[H,P g k] are deduced. 
The RPA equations for coordinate q q k and momentum p q k variables read 
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Here {ph\A q , k ,) is the matrix element for the two-quasiparticle state \ph), e p h is the energy of this 
state, ujy is the energy of the RPA state \v). The RPA phonon operator reads 



q pheq 



(11) 



where d h (d p h) are the creation (destruction) operators of two-quasiparticle configurations with 
the amplitudes determined by solutions of ©-([B]): 



pheq 



q'k' 



1{e v h ±v u ) 



(12) 



Following ©-(16]), the separable ansatz ([2]) explores the residual interaction of the Skyrme 
functional through the second functional derivatives. The calculations show that, for spin-flip 
magnetic modes, the spin 
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are most important. In deformed nuclei, the magnetic modes couple to electric motion 26)]. Then 
the relevant electric terms [lrj . \r\ should be added. For the scissors mode, which exists only in 
deformed nuclei, the spin-orbit, tensor, and electric terms are essential. 

As was mentioned above, the model is self-consistent in the sense that both the static mean 
field 



aq 



(17) 



and the residual interaction ([2]) are derived from the same functional. The rank of the RPA matrix 
©-© is determined by the number K of the input Q q k or P q k operators. Usually K = 2 5 



and so the rank is small [171 . l24l . |25| ] . This reduces the computational effort and allows systematic 
studies even for heavy deformed nuclei. 

The pairing functional reads V piL i r = l/2^2 q G q x q X q where Xq ls the pairing density and G q is 
the pairing strength [7]. In the present study, the pairing is included at the BCS level through the 
quasiparticle energies and Bogoliubov coefficients. 



B. Strength function 

GR in heavy nuclei are formed by many RPA states whose detailed structures cannot be resolved 
experimentally. Then a direct computation of the strength function is more efficient and reasonable. 
In SRPA the strength function for Ml excitations reads 

S{Ml-u) = Y,\^u\M)\ 2 Q{u-uj v ) (18) 
Bf-s + V \(ph\M)\ C(u-£ ph ) 

\z=u>H$ q % q 

where (^f u \M) is the matrix element of Ml transition between the ground and excited RPA 
states, and £(uj — u v ) = A/[2tt[(o; — u> u ) 2 + A 2 /4]] is the Lorentz weight to simulate the broadening 
effects beyond SRPA (escape widths, coupling with complex configurations). In the present study, 
the Lorentz averaging parameter is A=l MeV. Further, (3 = qkr with r labeling X and ^-operators, 
9 means the imaginary part of the value inside the brackets, F{z) is the determinant of the RPA 
matrix ©-([8]) with lu u replaced by the complex argument z, Fpp/{z) is the algebraic supplement 
of the determinant, and 

w avWjKp^M) n (Y) ( , ieph<ph\Yf k )(ph\M) 

U qk \ Z > ~ 2^ £ 2 _ 2 ' U qk \ Z > ~ 2^ £ 2 _ z 2 " 

q',pheq' P h q',ph€q' P h 



The operator of Ml transition reads M(Ml/i) = //s\/3/(87r) Ylq\9sSn + ff?^*] where spin ^-factors 
g§ = 5.58<; and g™ = —3.82? are quenched by ?=0.7 and orbital ^-factors are gf = 1 and gf = 0; 
and In are spin and orbital operators. 

In the present study, the strength function (|18|) is calculated only for the branch [i = 1 which 
gives the K n = 1 + states relevant for both spin-flip and orbital modes. 

C. Choice of input operators 

The SRPA formalism requires to choose for Ml modes the generating operators P q k(r). Their 
choice is crucial for a fast convergence of the separable expansion V^cF to the true Skyrme residual 



6 



interaction. We achieve this by using such P q k( r ) which allow the separable operators X q k(r) and 
Y q k(r) to have maxima in different spatial regions of the nucleus, in the surface and interior. Such 
a flexible choice allows to achieve good convergence already with a few separable terms. 

The physical arguments suggest that, for spin-flip mode, the leading input operator P q i(r) 
should have the form of the spin part of the applied external field (??) with \x = 1, i.e. P q \ = s q + . 
The detailed distributions depends on the interplay of surface and volume excitations. This can be 
resolved by taking into account the interior of the nucleus j uM by adding input operators 



with different radial parts, P q 2 = r s + , P q % = r s + . This results in X qk {r) and Y q k(r) interaction 
operators having more sensitivity in the interior. 

In deformed nuclei we should take into account the coupling between magnetic and electric 
= 1 + states. So the quadrupole input operator Q q = r 2 l2i with the counterpart P q = i[H, Q q ] 
should be added. This operator generates quadrupole excitations in both spherical and deformed 
nuclei. Besides it allows to explore the scissors mode in deformed nuclei. 

In Fig. 1, the spin- flip (g"' p = 0) strength function (|18p in 150 Nd is plotted for different sets of 
input operators. It is seen that all the sets give very similar results. In principle, it would suffice 
to use the minimal set P q = si and Q q = r 2 Y2\- However, to be on the safe side, we will use in the 
following calculations a set from 3 input operators, by adding also the operator P q = r 2 s q + . Then 
K=3 and we have the RPA matrix of the modest rank 4K=12. 




6 8 10 12 

CO [MeV] 

FIG. 1: Strength function for spin-flip Ml GR in 150 Nd, calculated with the force SkI3 for four sets of input 
operators, as indicated in the figure. 



Note, that in terms of two-quasiparticle matrix elements, the relations between input operators 
and their counterparts have the form 

Pqk -» Qqk(r) = i[H,P qk ] -» (ph\Q qk ) = 2e ph (ph\P qk ) - (ph\X q qk ) 
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for magnetic modes and 

Q q k -» P qk {r) = i[H,Q gk ] -> (p/»|^fc) = 2e pfc (p/i|Q, fc ) - (p/*^) 
for electric modes. For more details see fltl fis| . 

III. RESULTS AND DISCUSSION 

Results of the calculations are presented in Figs. 2-4. Fig. 2 shows the spin-flip and orbital Ml 
strengths for the isotopes 142_152 Nd. The strengths are computed with the force SkI3 by using gf=0 
and g s ' p =0, respectively. The isotopic chain ranges from semi-magic spherical (A=142) to axially 
deformed (A=150, 152) nuclei. The soft nuclei (A=144,146,148) require probably a beyond-RPA 
treatment. Nevertheless, they are included to illustrate the trends. 
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FIG. 2: Spin-flip (solid curve) and orbital (dotted curve) strength functions in 142_152 Nd, calculated with 
the force SkI3. For every isotope the parameter of quadrupole deformation (3 defined in Ref. |25[ is shown. 
For 152 Nd, the total strength (bold curve) with both spin-flip and orbital contributions is presented. 



Figure 2 includes the high-energy orbital strength at 30 MeV for all the isotopes. Actually 
this is the \fi = 21 branch of the isovector electric quadrupole GR, though with an overestimated 
energy. The low-energy scissor orbital mode is absent in spherical 142 Nd but then appears and 
steadily increases with growth of the deformation in heavier isotopes. In deformed 150 > 152 Nd, this 
mode dominates at 2-4 MeV and gives a significant contribution to Ml strength at 2-7 MeV. The 
comparison of the total, spin-flip and orbital strengths in 152 Nd shows that the interference between 
spin-flip and orbital modes can be both constructive and destructive. The results for other Skyrme 
forces used in the paper, SkT6, SkM*, and SLy6, are similar. 
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FIG. 3: Total spin-flip (bold curve), proton (solid curve), and neutron (dash curve) strength functions in 
142 Nd and 150 Nd, calculated with the force SkI3. The proton and neutron strengths are computed with 
<?" j = and g p s , = 0, respectively. 



Figure 2 also demonstrates a significant change of the structure of the spin-flip resonance with 
increasing deformation and neutron number. The distinct two-peak structure at 4 and 10 MeV in 
142 Nd evolves to a broad one-peak resonance in 150 > 152 Nd. Following our analysis of the structure of 



the RPA states, the spin-flip strength in Nd isotopes is determined by neutron u(lh^, 2 , 1/19/2) an d 



proton 7r(2cigy 2 , 2d 3 / 2 ) as well as ir(lg~^ 2 , I57/2) transitions. This is confirmed by Fig. 3 where pro- 
ton and neutron contributions are separated. One sees that in 142 Nd the neutron strength is given 
by one peak while the proton distribution exhibits two peaks, of 7r(2d7/ 2 , 2d 3/ / 2 ) and 7r(lc7^, 2 , L77/2) 
origin for 4 MeV and 10 MeV structures, respectively. The neutron contribution to the 10 MeV 
peak weakens from 142 Nd to 150 Nd since in the latter the neutron subshell I/2.9/2 is already partly 
occupied. Instead, the proton contribution is concentrated and enforced. 

Note that in Fig. 3 the proton and neutron peaks at 10 MeV are close to each other and 
actually form one broad resonance. This is a consequence of close neutron and proton spin-orbit 

23j. For the comparison, the forces SkM* 
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21 



22l | give a stronger 



splittings provided by SkI3 for medium and heavy nuclei 
and SLy6, with more different neutron and proton spin-orbit splittings 
splitting of the resonance in both 142 Nd and 150 Nd (Fig. 4b, d). Thus the Ml strength distribution 
is related to the spin-orbit structure of the underlying single-particle states. 

The above consideration indicates that the structure of spin-flip Ml mode depends on the 
Skyrme force. This is further illustrated in Fig. 4 where the results for four Skyrme forces are 
presented. The forces yield quite different results. Besides, none of these forces provides a proper 



271 ]. In particular, the high spike at 6 



quantitative description of the experimental data for 150 Nd 

MeV is not reproduced. This confirms the finding [l(3] that most of Skyrme parnruelcnza lions do 
not fully describe the spin-flip spectra. Nevertheless, our calculations allow to suggest a possible 
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origin of three Ml peaks observed in nuclides at the onset of rare-earth region [271 . |28( |. Following 
the above analysis, they can be considered as tt^cC^, 2d 3 / 2 )j n i^-9g/2' an d ^9/2) 

spin-flip transitions. 



> 





5 10 

CO [MeV] CO [MeV] 

FIG. 4: Total strength function in 142 Nd and 150 Nd for the Skyrme forces SkT6,SkI3, SkM*, and SLy6. In 
150 Nd the experimental data 2l| are depicted. 



IV. CONCLUSIONS 

The orbital and spin-flip Ml giant resonances in Nd isotopes with A=142, 144, 146, 148, 150, 
and 152 were investigated within the self-consistent separable random-phase- approximation model 
(SRPA) [if], . The calculations show the appearance of the scissors mode while moving from the 
spherical 142 Nd to deformed 152 Nd. The evolution of the spin-flip strength with deformation and 
neutron number was analyzed. The observed three-bump structure of this strength is explained 
by contributions of Tr(2d^ 2 , 2d 3 / 2 ), ^(lg^, 157/2), and i/(lh~^ 2 ,lh 9 / 2 ) spin-flip transitions. It is 
shown that the largest Ml peak in 142 Nd is mainly determined by the neutron contribution. The 
calculations reveal a still poor quantitative description of spin-flip Ml resonance by present-days 
Skyrme parameterizations. Further improvement of Skyrme forces is necessary, especially of their 
spin-orbit parts. 

Acknowledgments 

The work was supported by the DFG RE-322/12-1, Heisenberg-Landau (Germany - BLTP 
JINR), and Votruba - Blokhintsev (Czech Republic - BLTP JINR) grants. W.K. and P.-G.R. are 



10 



grateful for the BMBF support under contracts 06 DD 9052D and 06 ER 9063. Being a part of the 
research plan MSM 0021620859 (Education Ministry of Czech Republic) this work was also funded 
by Czech grant agency (grant No. 202/06/0363). P.V. is grateful for the FIDIPRO support. 



[1] M.N. Harakeh and A. van der Woude, Giant Resonances (Clarendon Press, Oxford, 2001). 

[2] Electric and magnetic giant resonances in nuclei, ed. J. Speth (World Scientific, Singapore, 1991). 

[3] F. Osterfeld, Rev. Mod. Phys. 64, 491 (1992). 

[4] N. Lo Iudice, Phys. Part. Nucl. 28, 556 (1997). 

[5] M. Bender, P.-H. Heenen, and P.-G. Rcinhard, Rev. Mod. Phys. 75, 121 (2003). 

[6] D. Vretenar, A.V. Afanasjev, G.A. Lalazissis, and P. Ring, Phys. Rep. 409, 101 (2005). 

[7] JR. Stone and P.-G. Reinhard, Prog. Part. Nucl. Phys. 58, 587 (2007). 

[8] T. Lcsinski ct al, Phys. Rev. C 76, 014312 (2007). 

[9] W. Zow et al, Phys. Rev. C77, 014314 (2008). 

[10] P. Vesely, J. Kvasil, V.O. Nesterenko, W. Klcinig, P.-G. Rcinhard, and V.Yu. Ponomarev, Phys. Rev. 

C80, 031302(R) (2009). 

[11] T.H.R. Skyrme, Phil. Mag. 1, 1043 (1956). 

[12] D. Vauthcrin and D.M. Brink, Phys. Rev. C5, 626 (1972). 

[13] Y.M. Engcl ct al, Nucl. Phys. A249, 215 (1975). 

[14] P. Sarriguren, E. Moya de Guerra, and R. Nojarov, Phys. Rev. C 54, 690 (1996). 

[15] RR. Hilton, W. Hdhcnbcrger, and P. Ring, Eur. Phys. J. Al, 257 (1998). 

[16] V.O. Nesterenko, J. Kvasil, and P.-G. Reinhard, Phys. Rev. C 66, 044307 (2002). 

[17] V.O. Nesterenko, W. Kleinig, J. Kvasil, P. Vesely, P.-G. Reinhard, and D.S. Dolci, Phys. Rev. C 74, 
064306 (2006). 

[18] P. Vesely, Ph. D. thesis, Charles University in Prague, Czech Rep., 2009. 

[19] N. Lo Iudice and F. Palumbo, Phys. Rev. Lett. 41, 1532 (1978). 

[20] F. Tondcur, M. Brack, M. Farine, and J.M. Pearson, Nucl. Phys. A420, 297 (1984). 

[21] J. Bartel et al, Nucl. Phys. A386, 79 (1982). 

[22] E. Chabanat et al, Nucl. Phys. A627, 710 (1997). 

[23] P.-G. Reinhard and H. Flocard, Nucl. Phys. A584, 467 (1995). 

[24] V.O. Nesterenko, W. Kleinig, J. Kvasil, P. Vesely, and P.-G. Rcinhard, Int. J. Mod. Phys. E 16, 624 

(2007); Int. J. Mod. Phys. E 17, 89, (2008). 

[25] W. Klcinig, V.O. Nesterenko, J. Kvasil, P.-G. Reinhard, and P. Vesely, Phys. Rev. C78, 044313 (2008). 

[26] J. Kvasil et al, Phys. Rev. C63, 054305 (2001). 

[27] H.L. Wortche, Ph.D. thesis, Technischen Hochschule Darmstadt, Germany, 1994. 

[28] D. Frekers et al, Phys. Lett. B244, 178 (1990). 



11 



